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We derive a special scalar field potential using the anisotropic Bianchi type I cosmological 
model from canonical quantum cosmology under determined conditions in the evolution to 
O^i anisotropic variables j3± . In the process, we obtain a family of potentials that has been intro- 






duced by hand in the literature to explain cosmological data. Considering supersymmetric 
quantum cosmology, this family is scanned, fixing the exponential potential as more viable 
in the inflation scenario V((/)) — Vq e^^'^'^. 
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'^. ■ I. INTRODUCTION 

K^ I One of the main problems of inflationary cosmology is to find a mechanism to derive in a natural 

j^ ■ way the appropriate scalar field potential in order to develop enough e-foldings of inflation. By 

natural, we understand a mechanism from which some theory provides a scalar field potential that 
offers the convenient features of inflation. In this work we derive a scalar field potential from 
supersymmetric quantum cosmology that gives these conditions. 

In a previous work, we determined scalar potentials from an exact solution to the Wheeler- 
DeWitt ( WDW) equation in the quantum cosmology scenario [Ij] , using as a toy model an homoge- 
nous and isotropic cosmological model. There we focus on solutions that may be relevant for the 
early universe constructed within of WKB approximation. Recently, these scalars potentials were 
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obtained using a local supersymmetric scheme [2]. Nowadays it is a common issue in cosmology to 
make use of scalar fields as the responsible agents of some of the most intriguing aspects of our 







universe |3l4lll| , such as inflation 
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131 ] , dark matter and dark energy [ij] . The natural derivation 
of a scalar potential is a challenge, posing the following question: What physical processes provide 
the adecuate scalar field potentials that govern the universe in determined epoch? To answer this 
question, we use the ideas of quantum cosmology to solve the Wheeler-DeWitt equation with a 
particular ansatz for the Bianchi type I universe wave-function. In this scheme, we obtain two 
possible scenarios, the first one with an scalar exponential potential V(0) = Vq e"*"^, and the second 
one giving a family of potentials, similar to those obtained in our previous work [ij]. It is interesting 
that in the first scenario the A parameter is not fixed by the quantum scheme, remaining as a free 
parameter of the theory. To fix it we invoke supersymmetric scale, using the tools of supersym- 
metric quantum cosmology in order to find most viable scalar potential for the inflationary epoch, 
in this scale. To do this, we applied supersymmetry as a square root of general relativity 15l4l8l|. 



in which the Grassmann variables are only auxiliary and can not identified as the supersymmetric 
partners of the cosmological bosonic variables. Therefore, we construct a family of scalar potentials 
treating the quantum solutions to anisotropic Bianchi type I cosmological model in the anisotropic 
variables /3+ and /3_ . The conditions we use give us a special structure for the scalar potential; By 
simplifying the Wheeler-DeWitt equation we obtain two cases: one in which both parameters I3± 
have hyperbolic trigonometric functions as solutions, and another where /5_ (/3+) have a trigono- 
metric (hyperbolic trigonometric) behavior. This potential is also a good candidate, depending 
on the parameter value, in order to study inflation, dark matter, dark energy or tachyon models 
[l9|. The transform Wheeler-DeWitt equation can be solved using a particular ansatz in the WKB 
approximation (Bohmian representation [2()]). This method has been used in the literature 21] to 
solve the cosmological Bianchi class A models, and in a particular, our result in the second case 
is similar to the one found in reference [l|] for the isotropic Friedmann-Robertson- Walker (FRW) 
cosmological model. On the other hand, the best candidates quantum solutions become these that 
have a damping behavior with respect to the scale factor, in sense that we obtain a good classical 
solution using the WKB approximation in any scenario in the evolution of our universe 
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231]. 



The supersymmetric scheme have the particularity that is very restrictive because there is more 
constraints equations applied to the wavefunction. So, in this work we found that exist a tendency 
for supersymmetric vacua to remain close to their semiclassical limits, because the exact solutions 
found are also the lowest-order WKB approximation, and not corresponds to the full quantum 
solutions found previously. 



II. THE WHEELER-DEWITT EQUATION 

On the Wheeler-DeWitt equation there are a lot of papers dealing with different problems, for 



example, Gibbons and Grishchuk 2j] asked the question of what a typical wave function for the 
universe is. In reference [25] appears one excellent recopilation of paper on quantum cosmology 
where the problem of how the universe emerged from Big Bang singularity can no longer be 
neglected. Also, an important approach to this problem is the wave function proposal in which 
the universe would be completely self-contained without any singularities and without any adges. 
Our goal in this paper deals with the problem to built the appropriate scalar potential in the 
inflationary scenario. 

We start by recalling the canonical formulation of the ADM formalism to the diagonal Bianchi 
Class A cosmological models. The metrics have the form 

ds2 = _(n2 _ NJNj)dt2 + e2^(*)e2ftjW Ju\ (1) 

where N and Ni are the lapse and shift functions respectively, il.{t) is a scalar and /3ij(t) a 3x3 
diagonal matrix, /Sy = diag(/3+ + \/3/3_, /3+ — \/3/3-, —2/3+), u^ are one-forms that characterize each 
cosmological Bianchi type model, and that obey dw' = |Cj]^a;J Aoj^, C-j^ the structure constants of 
the corresponding invariance group 2a]. The metric for the Bianchi type I, takes the form 

the total lagrangian density function is given by 

-C Total = -Cg + ^A + ^ma.ttev,(j) = V ~E. 0^ ~ 2A) + £matter,(;i; (3) 

we use as a first approximation a perfect fluid and a scalar field as the matter content, in a comoving 



frame 



26i, 



/^Totai = V^(R - 2A + 16^Gp + -g'^'^d^cPd.cl) + V(0)), (4) 



and using ([2]) we have 

J-Total — e 



6 ^^2 _ l^-^2 _ ^^2 _ ^^, ^ ^g^^^^ _ 2^^ ^ Vi^)N 



N N' ^ N' N"^ ' 2 

where we redefined the original scalar field as = -v/l2<^- 



(5) 



The corresponding momentas are calculated in the usual way 

+ dp+ N '^^ 12 + 

dC 12;S'_ oo ^- e-3f^^^^ 
n_ = — ^ = ^e^^ ^ /3- = Nn_ 



dp- N 12 

(9/: 12(/3 oo e-3f^ 

— — ^e — > w = 

N ^ 12 



n^ = — = -^6^^'^ = — T^Nn^ 



now writing ([5]) in canonical form £can = HqCi — N'H where % is the hamiltonian density function, 

Lean = n^ri + n+/3+ + n_/5L + n^v^ 

'n^ - n^ - n^ - n^ - e^^ [3847rGp - 48A + 12V{lp)\) 
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we obtain the corresponding Hamiltonian density function 



e 
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3f7 



U = ^^7^^ (n^ - n^ - n^ - n^ - e^^ [384^Gp - 48A + 12y(¥9)]) (6) 



when we include the energy-momentum tensor for a barotropic perfect fluid p = 7/9, we have 

n = — — (^n^ - n^ - tf - n^ + 48Ae6" - 3847rGM^e-3(T-i)^ - 12e^"V(v?)j . (7) 

Imposing the quantization condition and applying this hamiltonian to the wave function ^, we 
obtain the WDW equation for these models in the minisuperspace by the usual identification PqM 



by — i9qM in ([7]), with q'^ = (ri,/3+,/3_, 93), and following Hartle and Hawking 23] we consider a 
semi-general factor ordering which gives 



H^ 



^' +^ + ^ + il + Q^ + 48Ae6^ - 3847rGM^e-3(^-i)^ - l2e^^V{^) 



^ = 0, 



91^2 5^2 Qp2_ Q^2 ^QQ 

where Q measures the ambiguity in the factor ordering between the scalar function il. and its 
corresponding momenta. This equation is not easy to solve, first because we do not have the 
structure of scalar potential, and second, it depends strongly on the class of scenario we analyze 
with barotropic equation. In the follow, for simplicity we shall use the inflationary case, 7 = — 1 

Using the following ansatz for the wavefunction "^{QjCp, f3±) = e^°'^^+^°'^^~'B[Q,ip), we obtain 
a reduced WDW 

^2 d"^ ^ d 



xP^W^'^Wi^ "" <"'*'' - ='**'"°"- - ''^'^» + "" 



2 = 0, (8) 



where the constant c^ = af + a^ . 

Eqn ([5D can be written in compact form as 



aE + Q—-u{n,^,Xeffy^ 



0, 



(9) 



48A-3847rGM_i 



where the d'Alambertian in two dimensions is redefined as D = —5^ + 5^, X^fj 
is the effective cosmological constant, and the potential l]{Q,,ip,A) = e^^ [12Y{ip) — Xcs] — c^. 
To solve (IHl) we take the ansatz, which is similar to the one used in the Bohmian formalism into 

n 

quantum mechanics [20|] 



H(J^,<y9) = Ty(0,<y9)e-^(^''^) 



(10) 



where S{Q, (p) is the superpotential function. Eq ([9]) can be written as the following set of partial 
differential equations 



{VSf - u 



dS 



W{DS + Q— I +2VW-VS 



aw + Q 



dW 



0, 
0. 

0. 



(11a) 
(lib) 

(lie) 



where the first equation is the classical Hamilton-Jacobi equation, which plays an important role 
in this work. The different terms in this equation are 



vw-vs 



{dnW) (dnS) + {d^W) (d^S) , (V)^ ^ - {dnf + (d^)' 



Any exact solution complying with the set of equations (llla|llb]lllcp will also be an exact solution 
of the original WDW equation. Following reference [l(, first we shall choose to solve eqns (jllah 
and (jllbp . whose solutions at the end will have to fulfill with eqn (jllcp which play the role of a 
constraint equation. 

Taking the ansatz [27| [Sg] 

1 



S(0, if) = -e-'''g{(p) + c(bin + b2A(^), 



(12) 



7:2 5 +7:2 h£ -i2Vi^) + Xeff 






b2dg_ _ 
3 dif 



^30 



0. (13) 



with Aip = ip — ifQ, LfQ is a constant scalar field, bi arbitrary constants, Eq (jllap is transformed as 

2 

This equation is more dificult to solve, at this point we introduce the main idea of the paper to 
obtain the scalar potential family, which are strongly dependents to solutions for the anisotropic 
variables f3±. We include two steps to solve equation ([13]): 



1. First, consider that the second and third parenthesis are null, but maintaining that c 7^ 0. 



The first condition implies that bi = \J\ + bg, and the constants ai and a2 are real, the 
solutions for the anisotropic variables fi± can be considered as hyperbolic trigonometric 
functions. The second condition becomes an ordinary differential equation for the unknown 
function g((/7), yielding 

g(^) = goet^'^, (14) 

with go an integration constant and a = ^ = — ^ ^ with b2 7^ 0. The scalar potential 



function become, when we take the first parenthesis in eqn ()13p . 

V((/?) = (4A - 327rGM_i) + Voe'^^'^, (15) 

Qp.2 

with Vo = ^j^i° 2 • With these results, the superpotential function (fT2|) is 

S(f7, 99) = ^e^^et^^ + c (i^l + b^O + baAt^ J . (16) 

Sustituting (I16p into (lllbj) . the corresponding solutions for the function W in the form e^ ^'''^ 
become W = Exp [^(Q — 3)ri — ^aAc/?] , we develop the following wavefunction 



ai/3+ + a2/3- + -(Q - 3)J1 - -aAy? 



^ = Exp 
with the constraint on the parameter a < ±6. 



e ^ 



(17) 



In the literature 28|] the scalar potencial type V(i?i>) = e'^'^ gives a power law in the classical 
scale factor, considering the flat FRW cosmological model when A < — \/2. If we consider the 
extreme values for the a parameter (Q=0) and the corresponding transformation between 
99 —)•</), we obtain the special scalar potential V{4>) = Vq e ^^''', for standard inflationary 
model, this class of potential has the advantage that classical analytical solutions can be 
found and for appropriate values of the parameters, inflation can be obtained. 

2. In the second step, we consider that the constant c=0, implying that ai = =bia2, then the 
solutions for the wave function for the anisotropic variables (3± are considered trigonometric 
functions for the variable /3+ and hyperbolic trigonometric function for the variable /3_ . Thus 
the superpotential term (I12p has the simple form 

Sin,^) = -e^''g{ip), (18) 



and equation ()13p becomes an ordinary differential equation for the unknown function g{ip) 
in terms of the scalar potential V{(p, Xes) = V((/?) 

2 



Aeff 

12 ' 



dg 

dip 



9g2(v9) = 12/i2 V((^) 



12 



12/.2v(^^^^^)^ 



(19) 



this equation is similar to the one obtained in reference 



1[. It is not surprising that this 



equation is similar to eqn (12) in reference [l|, because the anisotropic Bianchi type I cosmo- 
logical model is the generalization of the flat FRW model. The last equation has several exact 
solutions, which can be generated in the following way. Consider that V = g^F{g), where 
F(g) is an arbitrary function of its argument. So, eq. (fT9|) can be written in cuadratures as 

1 /" ding 



Aif 



"2\/3 



! + /^'F(g) 



(20) 



In this way, we can solve the g{f) function, and then use the expression for the potential 
term V = g'^F{g) back again to find the corresponding scalar potential that leads to an exact 
solution to the Hamilton- Jacobi equation (jllap . Some examples are shown in Table I. 



TABLE I: Some exact solutions to eq. (|T9|) and their corresponding scalar potentials, where n is any real 
number and Vq is an arbitrary constant, different to step one. The third line is equivalent to this obtained 
by the first step. 



F(g) 





Vo 

F05-" {n ^ 2) 

hi g 

(ln.g)2 



1{V) 



Exp [±3Av3] 



^4^sinh(±3A^) 



62 



Ai^ 






2/n 






V{ip,Xee) 





Vo 



Vo exp (aAip) , a ^ ±4\/3J| + fiWo 



-ifPVoe-^^'f 



r^e^'' , r = i 



2\/3 



2(2-ri)/n 



r;=f 



ue2" , u = {^ly/ZAipf - ■^_ 



4^ 



'U(ip) 



2v3/iA(y 



In this way, the superpotential S(ri,(/9) is known. 
For solve (jllbp we assume 



^ ^ ^[z{Q.)+i^{^)] 



(21) 



we arrive to a set of ordinary differential equations for the functions z{Q) and uj{ip) 



^^-« = ^' 



z(!i) = ^l±in 



d^g dg doj 



+ 2-^— = 3{k + 3)g, 



dip"^ d(p dip 
Thus, the exphcit form for the function W becomes 

dip 



3/c f 



dip 



dJlng) 



(22) 



W = eyip^ — 



3 ./ d^{lng) 



2 i (8^9)' 



The constraint (jllcp can be written as 



9^w + (5<^a; 



,2 ^ 



Q^ 



(24) 



(25) 



and 



d^u 



3k 



3/i 



2 5<^[V((/3,Acflf)] 



25^(lng) ' (5^g, 

Taking into account table I, we present the corresponding wave function in each case, in table 



II. 



TABLE II: Wave function corresponding to table I 



g('^) 



wave function ^ 



3!!±3A.^ 



Exp [±3A^] 



exp {a2(±i/?+ + /?-) + (^) n ± f A(p} 



^i^sinh(±3A(p) 



cosli^(±3A(^) exp|a2(±i/3+ +^-) + (^) »} 



c='"V'^^4^..nh(±3A.^) 



efAv 



■27n 



exp {a2(±i/3+ + /?-) + (^) n + (^^34^) A^l 



c''^'^-4m''Voc~''^^ 
2\/3 



exp {a2(±i/3+ + /?-) + (^) ^i + w((^)} e" 



^3n+g(v5) 



a;(^) = I [-A,, + A Ln (V^Vo + e^''^.)] _ t^ Aretanh ^3^^^. 



,3ii+u(^) 



,"(¥>) 



exp 



{a2(±i/3+ + /3_)+('^)» + 



k+3-2u^ T A 3u^ A 2 

4^2 ^ LnA(y9 - -^Av3^ 



„3ii + r(y) 



■■(v) 



exp {a2(±i/3+ + /3_) + (^) fi + c^((^)} e" 
a;((^) = ^ Arctan (^e"^'^)] + 6^i^Aip - f Ln ( 



3 



„2u((p) 



_3 quI'P) _ ^-"(v) 



^Arctan (■^e""^'^: 



In special case, the third line in table II, have a damping term that correspond to e~^ and a 
plane wave type, similar to equation (J17p obtained in the first step. The first line has this behaviour, 
but corresponds to null scalar potential. 



In this way, using the quantum forniahsni in the sector of inflationary scenario, we found that the 
scalar potential becomes an exponential behaviour, however the coupling constant is undetermined. 
The question is, how can we fix the coupling constant?. The answer could be in the supersymmetric 
quantum cosmology using differential operators to the Grassmann variables, where we present the 
formalism in next section. 

III. SUPERSYMMETRIC QUANTUM MECHANICS 

In the following we shall apply the supersymmetric quantum formalism at the quantum structure 
obtained in the previous section, to obtain a closed value to the parameters that optimize the 
inflation scenario, i.e, we do an analysis to the family obtained in table I for the function g{(p); or 
in other words, which is the constraint on the superpotential function that appears in the quantum 
level, equation (llSp . based in tables I or II?. In this order of ideas, we found one integrability 
condition on the g{(p) function, which fix the coupling parameter of our problem. 

For obtain these results, in this section we consider only a reduced supersymetry in two bosonic 
variables {0,,(p), without consider the anisotropic parameters, due that the full problem do not 
contemplate the initial condition on our original problem. For instance the decomposition of the 
wave function in the full expansion have 16 components, or 16 x 16 matrix components, and 
the solution is very complicated. In this sense, for solve our problem we use reduced bosonic 
hamiltonian, eqn ^, and in consequence, one reduced supersymmetry. 

The idea of Witten [29( is to find the supersymmetric supercharges operators Q, Q that produce 
a superhamiltonial Tisusy, that satisfies the closed superalgebra 

^susy = ^ [Q, Q] , [Q, Q] = [Q, Q] = 0, (26) 

where the super hamiltonian Hsusy has the following form 



/i-susy — 'ib ~r 



^M 



r,v , (27) 



with %[) is the bosonic hamiltonian ([8]) taking the constant c=0, then the supersymmetric approach 
will only be applied to reduced hamiltonian, and S is the corresponding superpotential function 
that is related with the potential term that appear in the bosonic hamiltonian, i.e, have the same 
structure that in the quantum level, proposed in the last section. This idea was applied in reference 

n n 

[18l ] for all Bianchi Class A models without matter content, and in [30] to FRW cosmological model. 
For example in reference 31| it is explained that in particular, we can formulate a particle dynamics 
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in a potential V(q'^) on a curved manifold and super symmetry requires that the potential V(q^) 
is derivable from a globally defined superpotential S(q^) via V(q^) = |G^'^(q) ^y q^^ , where 
G'^'^(q^) is the metric in the curved space. This equation is represented in the quantum level by 
eq. (fTTi|) . 

In this approach, a supersymmetric state with Qlip >= is automatically a zero energy ground 
state, in a similar way that in the quantum regime. This simplifies the problem of finding supersym- 
metric ground state because the energy is known as priori and the factorization of T-Lsusyl"^ >= 
into Qlip >= 0, QIV' >= often provides a simpler first order equation for the ground state 
wavefunction. The simplicity of this factorization is related to the solubility of certain bosonic 
hamiltonians. In this work, as in others, we find for the empty (+) and filled (-) sector of the 
expantion of the wavefunction in this approach, in the sector of the fermion Fock space zero energy 
solutions \A± >= e^^|ib > where A± are the corresponding components for the empty and filled 
fermionic sector. 

The corresponding supercharges that satisfy the superalgebra, when we consider the bosonic 
hamiltonian given by equation ([8]) become 



Q 
Q 



r 



d dS 

+ 



d dS 



We consider the following algebra for the fermionic variables 



and the corresponding representation 



{v'^v'} = {V'^V''^} = o, 



(28) 



(29) 



w 



Equation (j28|) are 



^f- 



flU 



de^ 



Q = - 


' d dS' 

dqO ' aqO 


d 




' d as 

aql ' Oql 


' d 


Q = e° 


' d dS' 


+ e^ 


I 


d dS' 

?ql dq\ 





The decomposition of the wavefuntion becomes 

E{n, ip) = A+ + BqO^ + BiO^ + A^e^e\ 



(30) 



(31) 



(32) 



where the coordinates fields are q'^ 
contributions to the wavefunction. 



(q , q ) = (ri,c/?), A±, Bq,Bi are the bosonic and fermionic 
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The super symmetric equations Q|H >= 0, Q|S >= are 



QE 



+ 



d dS 
d dS ' 



d 



aO/iii 



^ [A+ + Bo9'' + Bi0' + A-e^e'] , 



Q: 



" d dS 

" d dS 
dq^ dq^ 

Then (j34p gives the following set of differential equations 



[A+ + BqO^ + Bie^ + A-e°0^] 



qO . 



'dA+ dS 

'dA+_ dS_ 

whose solutions to equations (I35I36P are 

A+ = a+ e^ 



0, 



0, 



aO/nl 



dq^ dq^ 



On the other hand, eq. ([55]) gives 
free term : 



^ + ^0— ' 
QqO QqO 

dA- dS_ 

'W^ dq^. 

dA- dS_ 

dq^ ^ dq^ 



+ 



dBi dS 

dq^ dq^ 



0, 
0, 



where (1391) can be written as 



(33) 

(34) 

(35) 
(36) 
(37) 

(38) 

(39) 
(40) 
(41) 

(42) 



V'"' {df,B, + B.d^S) = 0, 
considering the following ansatz for the fields B/j, 

B, = e^d,f+, (43) 

([37l) is satisfied identically, and (j42]l is 

rj^'" {d^dj+ + dj+d,,s + a,f+a^s) = r?^^ {d^dj+ + 2dj+d^s) = o, (44) 

where a possible solution is f+ = h.{il.zizip), and h is any function dependind to the argument, given 
the following constrain on the superpotential function 



dS__ dS_ 
dQ dip' 



(45) 
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and considering the structure of the superpotencial (jlSp we find one condition on integrabihty over 
the function g{^), given 

g((^) = goe±3A.,_ (4g) 

and taking into account tables I or II, we obtain the following constraint in the parameter a of the 
models when this last equation (j45p is satisfied 

- = ±3, 
2 

so, only the exponential scalar potential can survive in table I, and the coupling constant become 
a = ±6, given the scalar potencial V((?i>) = Vqc"^^ '^. In this way, supersymmetric quantum 
mechanics fix the values for the a parameter, being valid the argument introduced in the quantum 
scheme. 

Equations, (^0} HTj) can be written as 

with solution 

A- = a_e"^, (48) 

then, the set of contributions for the supersymmetric wave functions are find, 

A± = a±e^^ 
Bo = e^ao(/+) 
Bi = e'diiU) 

It is interesting to note that supersymmetry is very restrictive because exist more constraints 
equations applied to the wave function. In this sense, we observe a tendency for supersymmetric 
vacua to remain close to their semi-classical limits, because the exact solutions found are also the 
lowest-order WKB approximation. 

IV. CONCLUSIONS 

Using the quantum formalism in the inflationary scenario, we find that the scalar potential has 
an exponential behaviour as a good candidate. However, the coupling constant is undetermined. 
The question was, how can we fix the value of the coupling constant? The answer was in the 
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supersymmetric quantum cosmology using differential operators to the Grassmann variables, where 
the coupling constant is found under one condition of integrability on the function g{(p) = goe ''', 
and taking into account tables I or II, ^ = ±3. So, the main goal in this paper was to fix the 
value for the coupling constant to the inflationary scenario A = ^ = ±3 using the supersymmetric 
approach, when the quantum approach only gives the general structure for the scalar potential. 
Also we find exact solutions in both regimes. In the quantum level, we found that the possible 
solutions become the contributions to the empty(+) and filled (-) sector of decomposition to the 
wavefunction in the supersymmetric approach. 
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